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Abstract 

High sensitivity polarization measures, on wide angular scales, together with data 
on anisotropy, can be used to fix DE parameters. In this paper, first of all, we aim to 
determine the sensitivity needed to provide significant limits. Our analysis puts in 
evidence that there is a class of DE models that polarization measures can possibly 
exclude soon. This class includes models with DE due to a Ratra-Peebles (RP) 
potential. Using a likelihood analysis, we show that it is possible to distinguish RP 
models from ACDM and other dynamical DE models, already with the sensitivity of 
experiments like SPOrt or WMAP, thanks to their negative TE correlation at low- 
l, when the optical depth r is sufficiently large. On the contrary, fixing the energy 
scale A for RP potentials or distinguishing between ACDM and other DE potentials 
requires a much lower pixel noise, that no planned polarization experiment will 
achieve. While reviewing this paper after the referee report, the first-year WMAP 
data were released. WMAP finds large positive anisotropy-polarization correlations 
at low this apparently excludes DE models with RP potentials. 
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1 Introduction 



The nature of Dark Energy (DE) is one of the main puzzles of modern cos- 
mology. DE was first required by SNIa data (see, e.g., Perlmutter et al. 1998, 
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Riess et al. 1998), that indicate an accelerated cosmic expansion. A flat Uni- 
verse with fl m ~ 0.3 and VL^h 2 ~ 0.02 is also favored by an analysis of CMB 
and LSS observations (see, e.g., Tegmark et al., 2001; Netterfield et al, 2002; 
Percival et al., 2002; Efstathiou et al, 2002; Schuecker et al., 2003; Pogosian et 
al.,2003; here fl mt b- matter, baryon density parameters; h: Hubble parameter 
in units of 100 km/s/Mpc; CMB: cosmic microwave background; LSS: large 
scale structure). 

The cosmic component allowing flatness, in spite of Vl m ~ 0.3, should not be 
observable in the number-of-particle representation, in the present epoch. If 
this component is a false vacuum, whose pressure and energy density (pde 
and Pde) have ratio w = —1, a severe fine tuning of Pde, at the end of 
the electroweak transition, is required; in fact, before it, the vacuum energy 
density could be more than 10 50 times greater than now. A possible alternative, 
apparently reducing fine tuning requirements, is that DE arises from a self- 
interacting scalar field <fi (dynamical Dark Energy or quintessence; see, e.g., 
Ratra-Peebles 1988, RP hereafter; Wetterich 1988, 1995; Ferreira & Joyce 
1998; Perrotta & Baccigalupi 1999; Steinhardt et al. 1999; Zlatev et al. 1999; 
Albrecht & Skordis 2000; Amendola 2000; Brax, Martin & Riazuelo 2000). For 
dynamical DE, 

pde = ^ 2 +v{4>) , PDE = ^- 2 -v{4>) ) (l.i) 

V((f)) being the self-interaction potential, a is the scale factor and dots indicate 
derivatives in respect to the conformal time t. The ratio w becomes negative 
as soon as p k = (j) 2 /2a 2 < V . For p k jV ~ 1/2, it is w ~ —1/3 and dynamical 
DE approaches open CDM. Still smaller p k /V ratios allow to approach w = 
— 1 and a ACDM behavior. In order to work out w and its time variations, 
Friedman equations, in association with the equation of motion of 0, are to 
be integrated. 

Solutions clearly depend on the shape of V, which, in principle, is largely 
arbitrary. However, among potentials admitting a tracker solution, a particular 
relevance is kept by RP and SUGRA (Brax & Martin 1999, 2000) expressions 

V{4>) = A 4+a /<p a , V{<f) = (A 4+ 70 a ) exp(4vr0 2 /m;) , (1.2) 

which originate within the frame of Supersymmetric (SUSY) theories. In the 
expressions (1.2), m p is the Planck mass and A is an energy scale, currently set 
in the interval 10 2 -10 10 GeV, ranging between the electroweak scale and the 
SUSY soft breaking scale. The potentials (1.2) also depend on the exponent 
a; once A and Qde are assigned, however, a is univocally determined. The 
RP and SUGRA potentials are also indicative of the opposite cases, when w 
shows slow or fast time variations, respectively. 
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Dynamical DE and ACDM predict quite similar observational results, in most 
cases. This is a point in favor of dynamical DE, as it is known that ACDM 
allows a good fit of most available data. It is however important to devise 
suitable experiments whose outputs enable us to probe DE nature; if DE is 
dynamical, they should also enable us to fix the expression of V and the values 
of the parameters on which V depends. 

Here we show that measures of anisotropy and polarization of CMB at large 
angular scales can constrain V. If favorable (but not exceptional) conditions 
are fulfilled, polarization measures, at a sensitivity level comparable with ex- 
periments in progress, can provide significant information, possibly excluding 
that a RP potential describes DE. 

We shall show this point by performing a likelihood analysis, assuming that 
polarization measures come from an experiment similar to the Sky Polarization 
Observatory (SPOrt: see, e.g., Macculi et al, 2000; Carretti et al, 2000; Peverini 
et al, 2001), considering both a noise level that the experiment will attain and 
a lower noise level, as well. Our results, however, are not directly related to 
the above experiment and can be easily extrapolated to other observational 
contexts. 



2 The cosmic opacity to CMB photons 

In the present cosmic epoch, although diffuse baryonic materials are almost 
completely ionized, the scattering time for CMB photons t s ~ 4.45-10 18 r2^" 1 /i _2 s 
exceeds the Hubble time tu = 3.09 • 10 17 h~ 1 s and the Universe is transparent 
to CMB photons. These figures can be extrapolated to the past, but, while 
t s oc a 3 , the dependence of tu on a depends on the cosmological model. In 
an expansion regime dominated by non-relativistic matter, a fully ionized 
Universe is opaque to CMB at redshifts z > z op ~ 5(Q b hy 2/3 . Therefore, 
assuming Vl^h 2 ~ 0.022 (standard BBNS limits; see, e.g., Dolgov 2002 for a 
recent review) and h ~ 0.5-0.7, we have z op > 64, safely above any expected 
value for the reionization redshift z r i. No reasonable cosmological model allows 
a substantial reduction of such z op . 

While, therefore, r <C 1, its specific value depends on the model and, in 
particular, on the reionization history. If reionization is a single and sudden 
event, it can be synthesized by a value of z T {. In Fig. 1 we show how r is related 
to z ri , for a wide set of cosmological models. In general, z ri ~ 8 or ~ 16, if we 
require r = 0.08 or 0.16. 

Data can constrain the reionization history in indirect ways. High-z QSO 
spectra (Djorgovski et al. 2001; Becker et al. 2001) show that at least a fraction 
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of neutral hydrogen is present, in the intergalactic medium, at z > 6. These 
very authors, however, warn the reader against an immediate conclusion that 
z ri ~ 6; an ionization rate x ~ 0.01 can cause the Gunn-Peterson effect 
observed. Their finding, however, was taken as an indication of the expected 
range of z T {. 

Another constraint can be found, if the ionizing flux observed today is ex- 
trapolated to high z; in this way, Miralda-Escude (2002) found that that the 
photon flux is able to reionize the cosmic hydrogen at z ri ~ 8-10 and r ~ 0.08. 

Pushing farther the theoretical analysis, Tegmark (1997) estimated that, at 
z ~ 30, a fraction ~ 10~ 3 of baryonic matter is in fluctuations of ~ 10 2 -10 3 M Q , 
already in non-linear growth. They can cool and turn into massive stars thanks 
to molecular hydrogen radiation. However, as soon as thermonuclear processes 
ignite, stars produce photons hard enough to turn all hydrogen into atomic 
(or plasma) state and star formation is stopped until much greater fluctua- 
tions (~ 10 7 M Q ) enter a non-linear regime (see, e.g., Haiman 1997, Gnedin 
& Ostriker 1997). However, Haiman et al. (2000), Oh (2001) and Venkatesan 
et al. (2001) find that the situation can be even more complicated, as X-ray 
heating might enhance the formation of molecular hydrogen. From such pic- 
tures, Cen (2002) and Whyithe & Loeb (2002) try to deduce the reionization 
history. According to them, when the first PopIII stars form, a first high-z 
reionization occurs. PopIII stars eventually produce metals that, in turn, allow 
less massive stars to form and, therefore, the ionizing photon flux is reduced 
soon. Then the overall ionization rate is set by a balance between star for- 
mation and hydrogen recombination; the former overwhelming the latter at 
a redshift ~ 8, when x safely approaches unity. At such redshift, this picture 
meets Miralda-Escude (2002) conclusions and is therefore also coherent with 
the ionizing flux observed today. 

These pictures are qualitatively realistic, although their quantitative aspects 
are still uncertain. It should be emphasized that different reionization histo- 
ries, also when they yield the same r, give different CMB angular spectra. In 
this work we shall quantitatively explore simple reionization histories, where 
a single reionization event yields r = 0.08 or 0.16 and comment on the exten- 
sion of our results to more complex histories. The former r value agrees with 
Miralda-Escude (2002) estimate. The latter value assumes an earlier reioniza- 
tion, implicitly allowing a small hydrogen fraction (~ 1 %) to remain neutral 
until z ~ 6, to account for the observed Gunn-Peterson effect in high-z QSO 
spectra. 

It should also be reminded that, until CMB spectra analysis is restricted to 
anisotropy data, no significant constraint on r can be set; Stompor et al (2002), 
using balloon data, can only set the limit r < 0.40, at the 2-a level. Limits 
are so loose because of the degeneracy between r and n (primeval spectral 
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index), in CMB anisotropy data. In a recent paper, Colombo & Bonometto 
(2003) discussed how much such degeneracy can be removed using polarization 
measures. Recent observations (Kovac et al. 2002) apparently confirmed that 
CMB is polarized, but the sensitivity level and the angular scales probed by 
such observations were not adequate to reduce the above n-r degeneracy. 

In this work we explore what large angular scale CMB polarization data can 
tell us on DE nature. For the sake of simplicity, we assume that all model 
parameters are obtained from measures at smaller angular scales and concen- 
trate on the effect of the potential V(<f>) on large angular scales, when r is 
large. 



3 CMB angular spectra 



The linear evolution of density fluctuations in models with a DE component 
yields the angular spectra Cf; A = T, E, B, X stand for anisotropy, E- and 
-B-mode polarization and anisotropy-polarization cross correlation (X=TE) 
only E-mode polarization is considered through this paper. In Fig. 2-5 we 
show the spectra C l ' ' when the potential is RP or SUGRA and r is 0.08 or 
0.16, for a set of values of the energy scale A. All spectra are obtained with 
suitable extensions of CMBFAST 4 . 

Figs. 2 to 5 show that the nature of DE and, in particular, the value of the 
energy scale A, both for RP and SUGRA models, display their main effects 
on low-/ TE correlation spectra. In particular, Figs. 2 and 3 show that is 
negative, at low /, for RP potentials, while, in SUGRA models, such feature is 
absent. Negative are not an exceptional feature: all models yield alternate 
intervals of positive and negative C ; x as I increases. In most spatially flat mod- 
els, however, the spectrum, at low / values, starts from positive values. On 
the contrary, open CDM models, when r is large enough, start with negative 
Cf ', al low I. 

In the discussion section, we shall discuss why RP models share this feature 
with open models. This arises from the way how p DE scales with a, similar to 
the scaling of the curvature term in open models. On the contrary, in ACDM 
models, p^E does not depend on a. When the scaling of pde with a is mild 
and, therefore, close to ACDM, no low-/ negative Cf is expected. 



4 http : //physics . nyu . edu/matiasz/CMBFAST/cmbf ast . html 
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4 Likelihood analysis 



4-1 The SPOrt experiment 

In this section we aim to determine the likelihood of different models in re- 
spect to possible results of anisotropy and polarization measurements, ob- 
tained through definite observational apparati, for specific cosmological mod- 
els. In particular, polarization data sets are built with reference to the features 
of the SPOrt experiment. 

At variance from other space experiments (WMAP and PLANCK), SPOrt 
has been explicitly designed to measure the Stokes parameters Q and U and 
to minimize systematics and instrumental polarization. SPOrt will observe a 
sky area with declination \8\ < 51.6° (~ 80% of the whole sky). The lifetime 
of the experiment will be at least 1.5 years, but an extension is likely. 

Here we assume that pixels are distributed according to the HEALPix 5 pack- 
age, with N slde = 16, and smoothed with a Gaussian beam of FWHM of 7°. 
We therefore have, on the whole sky, 3072 pixels, whose centers lie at an av- 
erage angular distance of ~ 3.7°. Once polar caps are excluded, there remain 
2448 pixels, providing measures of the Stokes parameters Q and U. 

SPOrt polarimeters will provide no anisotropy data. Such information will 
be however available through other experiments (COBE, WMAP). Of course, 
anisotropy data probe a different sky area, i.e. the whole sky, excluding the 
area with galactic declination |<5| < 20°, where contamination by the Milky 
Way is severe. (On the contrary, such contamination can be assumed to be 
under control, for polarization data; see, e.g., Bruscoli et al. 2002; Tucci et al. 
2002.) In this sky area there are 1984 pixels. For 1360 of them, both anisotropy 
and polarization data are supposed to be available. 

Random noise (<Jp ix ) is included in artificial data, assuming it to be uncor- 
rected, both among pixels and between Q and U, in the case of polariza- 
tion measures. We assume a^ ix = 1 //K for temperature data (as expected for 
WMAP measurements, extended to 4 years, scaled to the SPOrt resolution). 
On the contrary, we consider different noise levels for polarization data. 

At present, the SPOrt experiment team expects a pixel noise level of 



where A is the experiment's lifetime in years while e indicates the detection 
efficiency Our analysis considers crf^'s in the interval comprised from 0.1 to 




5 http : //www . eso . org/science/healpix/ 
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1 /iK per pixel. Even at its highest value, this noise is consistent with SPOrt 
claims only if A = 4y and efficiency approaches unity. However, the SPOrt 
payload is to be delivered in more than one year and there is room for some 
further improvement in the apparati. 

4-2 Models explored 

In this analysis we take Vl DE = 0.7, VL b h 2 = 0.022, h = 0.7 and a primeval 
spectral index n—1. We considered also different values for these parameters, 
but the likelihood distributions are essentially independent from them. On the 
contrary, such distributions depend on the choice of r; as already outlined, we 
studied likelihood distributions for r = 0.08 and 0.16 . 

We consider first RP models, for A/GeV= 10 3 and 10 6 . As is known, RP 
models are accelerated only if 10 4 GeV and disagree with SNIa data at the 
2-er level, when A> 10 6 GeV. On the contrary, when A^ 1 GeV, the angular 
spectra of the model are hardly distinguishable from ACDM. Accordingly, 
fits yielding A^ 1 GeV, for a RP model, indicate that data do not allow to 
appreciate that the model is different from ACDM. 

We consider then SUGRA models, for A/GeV= 10 3 , 10 6 and 10 9 . Their TE 
correlation spectra show differences from ACDM and among themselves, al- 
though not so relevant as for RP. 

TE correlation spectra also depend on the reionization history. In this paper, 
such dependence is not explored and this is the main quantitative limitation 
to the results shown here. 

Dealing with low Z's, cosmic variance must be carefully taken into account. We 
explore this point by simulating and analyzing a large number of artificial data 
sets. Results will be both synthesized into average predictions and analyzed 
providing the frequency of possible determinations. 

4-3 Results 

We report first results for RP models, for A/GeV= 10 3 and 10 6 . In Figs. 6 
and 7 we show the average likelihood distributions for such models. In Figs. 8 
and 9, we report the distributions of the 1-2 a lower limits on A, for r = 0.08 
and 0.16, A = 10 3 and 10 6 , a? x = 0.1 //K and 1/iK. 

Fig. 10 is analogous to Fig. 6, but concerns a SUGRA model with A = 10 9 GeV 
and r = 0.16; among the cases considered here, this is the one which allows 



7 



the best recovery of A. In Fig. 11, we report the distributions of the lower 
limits on A for such model (at 1 a) . 

The above figures show that the energy scale A cannot be detected for SUGRA 
models, at least with large-angular-scale experiments. This is not only true at 
a noise level achievable by SPOrt, but also when noise is reduced by a factor 
~ 10. 

The situation was better for RP models. The capacity of an experiment to 
distinguish a RP model from ACDM increases with the energy scale A and 
the opacity r, although being more sensitive to the former parameter. In fact, 
a substantial fraction of realizations, with A = 10 6 GeV and r = 0.08, can be 
distinguished from ACDM, even at the 3-cr level; such fraction is only slightly 
greater for r = 0.16. Furthermore, only in quite a small fraction of realizations, 
RP cannot be distinguished from ACDM, at the 2-a level. Both with r = 0.08 
and 0.16, also models with A = 10 3 GeV can be distinguished from ACDM in a 
smaller (but still substantial) fraction of realizations, at least at the 2-cr level. 

If we aim farther and try to determine A, we meet substantial difficulties: even 
reducing the noise level by a factor ~ 10, the best we can find is a lower limit 
on the order of magnitude of A. Unfortunately, this can be fixed just in a 
fraction of realizations. 



5 Discussion 

Let us first discuss how the likelihood distribution can be obtained from the 
spectra of a given model, taking into account that the number of pixels for 
anisotropy and polarization (N T and N P ) in our (artificial) data are however 
different. Let then Tj be the anisotropy measured in N T pixels and Qj and Uj 
the Stokes parameters measured in Np pixels. In general, let us define vectors 

x = (Ti, , T Nt , Qi, , Q Np , Ui, , U Np ), of N s = N T + 2N P components, 

defining an observed state of anisotropy and polarization. Once a A value is 
assigned, the angular spectra Cf 4 = (Cf, Cf , Cf) are univocally determined. 
On the contrary, a data vector d, of N s components, built from them, is just 
a realization of such model: once the N s component vector d is assigned, the 
value of A is not univocally fixed. 

A function 

L(d\Cf) oc [detM]-^exp[- ^d T M-M] (5.1) 

is then built, to yield the likelihood of a given set of Cf (i.e., of a given A 
value), if d is observed. The main ingredient of L is the correlation matrix 
My = (xfxj) = Sij + Nij] here S^- is the signal term and is due to 
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the noise. The components Mij yield the correlation between the ith and jth 
elements of data vectors x corresponding to particular choices of C/ 4 , i.e. of A 
values. The construction of the (model dependent) signal term, however, does 
not require to build explicitly the vectors x. The procedure to be followed, in 
the case when both anisotropy and polarization data are available, is explicitly 
reported by Zaldarriaga (1998). The construction of the noise term is simpler, 
as we expect no noise correlation, and the matrix = 5ij &T,pix 2 (for % = 
1, N T ) and N i3 - = 8^ <Jp tPix 2 (for i = N T + 1, N s ) is diagonal. 

In what follows, the technical role of the Cf spectra does not need to be 
further outlined and the likelihood function is explicitly considered to depend 
on A. In particular, we assume that the most probable A value, for a given 
anisotropy-polarization state d, is the one which maximizes the likelihood. 
By integrating the likelihood distribution along the A axis, we are then able 
to find the intervals corresponding to 1-2-3 a levels. 

In order to explore cosmic variance, 1000 realizations of each model were 
considered; the likelihood curves plotted in the Figs. 5-6-7-11 are averages of 
the results of such realizations. 1-2-3 er value distributions, instead, allow to 
appreciate the spread among the results of each such realization. 

The feature which explains why RP model likelihoods can enable us to dis- 
tinguish them from ACDM, in several cases, is the negative TE correlation at 
low I. In the case of SUGRA models such TE anticorrelation is absent and 
large angle data seem unsuitable to provide information on DE nature. 

Negative TE correlations are essentially related to the simultaneous action 
of ISW effect and cosmic opacity. The former effect is a consequence of the 
varying rate of expansion, when we pass from matter to curvature dominance, 
in open models, or from matter to vacuum dominance in ACDM models. 
However, while the former kind of passage, in the presence of cosmic opacity, 
does produce anticorrelation, the latter kind of passage does not. 

The capacity of RP models to induce anticorrelation, as we shall see, is related 
to their features, which more closely approach open CDM, rather than ACDM. 
The opposite is true for SUGRA models. 

A full understanding of model dynamics requires that the whole set of Boltz- 
mann equations, describing photon distribution, is followed in time. Such 
equations are shown in a number of papers, together with equations ruling 
the dynamics of fluctuations in other components of a model (see, e.g., Ma 
& Bertschinger 1995). Signs and other definitions that we adopt here are the 
same as in the public code CMBFAST. Let us then indicate by Fi(k,t) and 
Gi(k,t) the Boltzmann components for anisotropy and polarization, respec- 
tively, k and t being the wave-number and the conformal time. 
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Let then n e and oy be the free electron density and the Thomson cross- 
section. The opacity to CMB photons reads then r = J* t *° n e (t')aTa(t') dt' and 
— f = an e OT- 

Let us report then the set of equations fulfilled by the G\ components in 
spatially flat models: 

Gi = -f [ -Gi + ^(5 l0 + S f)} + JL^[l G,_x - (I + l)G l+1 ] (5.2/) 

(Sin is the Kronecker symbol). Here 

IT = G Q + G 2 + F 2 , (5.3) 
is the only vehicle transferring signals from anisotropy to polarization. 
In a non-flat geometries eq. (5.2f) becomes 

G l = -r [ - G/ + 5(5/o + y)] + 2^iP 6 ' G '-i " ( Z + ^^i^+i] J ( 5 - 2 °) 

here (3 2 = k 2 + K and bf = 1 - fT/ 2 //? 2 , with X = -(1 - n m )H 2 , H being 
the present value of the Hubble parameter. 

Quite in general, when n e and, therefore, — f is large, all Gi are exponentially 
damped. On the contrary, when n e is low, the polarization signal spreads along 
the Gi component series. 

In order to produce polarization, therefore, there must however be a seed 
coming from anisotropy, linked to the non-vanishing of the quadrupole term 
F 2 (k,t). If F 2 (k,t) becomes significant only when n e is very low, little polar- 
ization is produced on harmonics with I corresponding to such k. It is so for 
wavelengths 2n/k entering the horizon well after recombination, which corre- 
spond to spectral components with I <ti 200, unless reionization occurs. This 
explains the sudden rise of Cf in Fig. 2-5, when r > is considered. 

Let F° (Gf) be the value taken by Fi (Gi) at the present time t a . In spatially 
flat models, the angular spectra read 

Cj = -J d*kP (k) |F°(A:)| 2 , C? = -J d 3 kP Q (k) \G°(k)\ 2 , 

C ? = \\ d 3 kP (k)F l °(k)G^(k) , (5.4/) 

P (k) being the primeval fluctuation spectrum. In open model, angular spectra 
expressions are slightly more complex and read 

Cf = -J <PpP (q) IW)| 2 , Cf = -J d^P (q) \Gm\ 2 , 
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Cf = \\ d 3 pP (q) F?(0)G>M ; (5.4o) 
here q = (f3 2 - 4K) 2 /[3([3 2 - K). Clearly, for Vt rn = 1, both /3 and q return fc. 

Both eqs. (5.2) and (5.4) indicate the presence of suitable shifts in the /c-space, 
when open models are considered. In particular, the b\ coefficients in eq. (5.2o) 
are responsible for the shifts of C\ peaks, while the shifts due to the passage 
from k to j3 and q, in eq. (5.4o), displace the power of the spectrum P Q through 
the harmonics Fi at small /, i.e. on scales comparable with the curvature scale. 

Apart of these geometric and power shifts, the main difference between flat 
and open models resides in the equation fulfilled by the gravitational field 
fluctuations h and r\. In the synchronous gauge, for open models: 

2k 2 f) = SnGa 2 [£(Pc + Pc )9 c - h Po , cr {l - Q m )/a 2 }; (5.5o) 

c 

here p 0jCr being the critical density at z = and J2 C indicates a sum over all 
(relativistic or non-relativistic) matter components. In the same gauge, for 
flat models with DE: 

2k 2 f) = 87rGa 2 [J2(Pc + Pc)O c + {p DE + Pde)O de }, (5.5/) 

c 

where the DE term has been set in evidence, with its possible fluctuations. 
The wave-numbers k 2 and k 2 are shifted by 3K. Once again this corresponds 
to a shift in the /c-space. 

Let us notice soon that, in a ACDM model, no DE fluctuations exist and, 
furthermore, p^E = ~Pde, so that the second term at the r.h.s. of eq. (5.5f) 
vanishes. 

In models with dynamical DE, instead, 6 DE does not vanish and, in general, 
Pde 7^ —Pde- Accordingly, the second term in square brackets in eq. (5.5f) may 
read 0DEPo,cr{^ — ^m){^ + w)(pDE/ Po,de) (the ratio in the last parenthesis tells 
us how DE energy scales with a). This term is then analogous to the second 
term in square bracket in eq. (5.5o). It would coincide with it if w = —1/3 
and, namely, if 9 = —h/2. If the latter is true, we can expect that, apart of a 
different power distribution along the / axis and geometric effects at greater 
I, there can be specific similarities in the behavior of open and dynamical DE 
models. 

The relation between h/2 and Ode can be studied through the equation 

Ode + \ = [Sde + S-(c 2 s - w)5 DE ], (5.6) 

whose validity indicates that DE behaves as a fluid; here c s is the sound 
velocity in DE. This equation can be derived from the equation of motion of 
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0, as shown in Appendix A. An order of magnitude estimate, however, tells 
us soon that Ode ~ k 2 tSr>E', then, the ratio between Ode and the r.h.s. is ~ 
(t/L) 2 , where L is the scale related to k. When t <^ L (before horizon crossing), 
the D e term is negligible, in comparison with the r.h.s.. For L <C t (after 
horizon crossing), the contrary is true. Hence, before horizon crossing, h/2 is 
essentially required to equate the r.h.s. and, therefore, the ratio —h/20DE shall 
exceed unity. At horizon crossing such ratio must approach unity and keep to 
such value as t grows greater than L. 

This point can also be numerically inspected. In Fig. 12, the ratio —h/20DE 
is plotted against a, for those k values which yield the top contribution to the 
/ harmonics indicated in the figure, in the case of a RP model. Clearly, about 
horizon crossing, DE is already quite close to —h/2. 

The main differences between open and dynamical DE models, in the r.h.s. of 
eqs. (5.5) are therefore relegated to times before horizon crossing. Afterwards, 
the residual difference is due to a factor 1 +w. If w is too close to — 1, however, 
such factor risks to spoil the similarities noticed hereabove. 

The changes in rj directly act on F 2 . The equation fulfilled by this spectral 
component reads: 

+ j[2F 1 (k,t)-3F 3 {k,t)] + 


+ |(«M) + ^) (57/, 

in flat models, while in open models we have: 



F 2 (k, t) = an e OT 



-F 2 (k,t) + 



U(k,t) 
10 



F 2 ((3,t) = an e cr T 



-F 2 ((3,t) + 
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^[2b 2 F 1 ((3,t)-3b 3 F 3 (P,t)}+ 



8(3 



b 2 [f ) (f3,t) + 



h(p,ty 



(5.7o) 



5 k * V" ' ' ' 6 

Let us recall that F 2 is the only term of II where anisotropy contributes. Such 
contribution is therefore directly affected by a change in the equation of the 
gravitational field, which describes the ISW effect in the synchronous gauge 
and can be similar in open models and in models with dynamical DE, provided 
that w is not too far from -1/3. 

Accordingly, when open models show negative TE correlations at low /, we 
can expect the same to occur in flat RP models. On the contrary, SUGRA 
models, which grant a cosmic acceleration much closer to ACDM models, due 
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to a w value not too far from -1, should not be expected to give negative TE 
correlation and cannot be tested in the same way as RP models. 



6 Conclusions 

In this paper we showed that a particular class of dynamical DE models can 
be fairly safely falsified using large angle CMB polarization data. Models with 
DE due to a scalar field, self-interacting through a RP potential, belong to 
such class. In general, however, we expect that any model with dynamical DE, 
yielding w ~ — 1/3 or just slightly smaller, can be tested through similar data. 

Our treatment provides quantitative predictions for the case of RP potentials, 
with the simplest possible reionization history yielding r = 0.08 or 0.16. If the 
reionization history was more complicated, as is however likely, quantitative 
conclusions do not apply and is hard to parametrize the whole functional 
parameter space, providing precise predictions. An analysis of TE correlation 
spectra, however, shows that negative values are a generic feature, at low I, 
for most reionization histories. Accordingly, in the absence of such negative 
signal, a RP model is however disfavored. 

An experiment to detect CMB polarization with a resolution of ~ 7° and 
a pix — 1 i s f un y realistic. With such noise, there is a substantial fraction 
of realizations of a RP model with A = 10 6 GeV, with r = 0.16, which can be 
distinguished from ACDM at the 3-cr level. Accordingly, the data analysis of 
a polarization experiment, with a noise level of ~ a^ ix = 1 fiK, has, at least, 
~ 30 % of probability to indicate that the model is RP, if this is true. At the 
2-cr level, such model could be confused with ACDM only in quite a small 
fraction of realizations. 

The only planned polarization experiment specifically aiming to large angle 
polarization data is SPOrt. The noise level that SPOrt now expects is not 
much more than 1 fiK per pixel. Other experiments, planned to provide data 
with greater resolution, can also provide large-angle harmonics. For instance, 
the WMAP satellite, although not built to minimize systematic effects on 
polarization, can be able to provide polarization data with a noise level just 
slightly worse than SPOrt, so allowing an earlier test of RP models. Finally, 
with the sensitivity expected for PLANCK, RP and similar models can be 
surely tested. 

In order to determine the energy scale A for RP potentials or to distinguish 
between ACDM from other DE potentials, using large angular scale data, a 
pixel noise even below 0.1 /iK is apparently required. No planned polarization 
experiment will achieve such noise level soon and we can expect that an in- 
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sight on the shape of DE potentials can come earlier from different kinds of 
experiments. 



7 Note added after WMAP release 

After submission of this paper and during its revision following the referee 
report, the first release of WMAP data occurred 6 . This suggests to add some 
further comments. 

Apparently WMAP data give no indication of being negative. The three 
procedures, used to estimate it, all yield large Cf values at low I and r = 
0.16 ± 0.04, for ACDM models and a single reionization. One of the above 
procedures explicitly assumes that a fair approximation to the true model 
can be found among ACDM models, trying to fix ACDM parameters best 
fitting data. The Cf spectrum was also related its Fourier transform, and then 
determined with a third procedure claimed to be fully model-independent. 

Here we wish to outline that the first procedure is unsuited to exclude RP 
models. Following the pattern of Sec. 5, we created a set of artificial skies 
for RP models with A = 10 3 GeV and 10 6 GeV and r = 0.16. We applied the 
maximum likelihood procedure to these data and fit cosmological parameters, 
assuming that the true model is ACDM. For the sake of example, in Fig. 13 
we show the likelihood distribution on for a RP model with A = 10 3 GeV. 
The procedure, being unable to explore a fair class of models, systematically 
overestimates floE] if the noise is ~ 1 /xK, the likelihood distribution is sig- 
nificantly non-Gaussian. When noise is reduced, the likelihood distribution 
is more Gaussian, but still peaks at high Vl DE . Then, for a^ ix = 0.1 /xK, the 
Gaussian becomes very narrow and the fair Qde value lies at more than 3-cr 
from the best fitting Qde- The situation is similar for A = 10 6 GeV. 

In Fig. 14, we make use of this latter model to compare the RP spectrum 
used to create data, with the best-fit ACDM spectrum. The continuous line 
connects the Cf values of the RP model; the long-dashed line yields the Cf 
values for a ACDM model with the right Qde', the short-dashed connects the 
values for the best-fit ACDM model, with \—a error-bars. The shaded 
area shows the range of the cosmic variance. 

For the sake of comparison, Figs. 15 and 16 are analogous to Figs. 13 and 14, 
but for a SUGRA model. In the absence of the peculiar negative behavior, 
characterizing the RP case, no misleading effect is present. 



http : //lambda . gsf c . nasa . gov/ product/map/ m_overview . html 
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With the reserves that these tests suggest, and relying on the model indepen- 
dent deduction procedure, we however conclude that WMAP apparently 
excludes that DE is due to a scalar field, self-interacting through a RP po- 
tential with A> 1-10 GeV. We plan to deepen the quantitative aspects of this 
finding in a forthcoming paper. 
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Appendix A 



Starting from the Lagrangian density for the field, self-interacting through 
a potential V(0), and using Eulero-Lagrange equations, we can work out the 
equations for the background unperturbed field and its fluctuations 5(f). The 
latter reads 

50 + 2-50 + (k 2 + V"a 2 ) 50 + 0^ = (Al) 
a 2 

and we aim to show that this is equivalent to eq. (5.6), that we report here: 

e + ! = -IT^ + 3 ^ • {A2) 

An equation equivalent to eq. (A2) is given, for a generic fluid, by Ma & 
Bertschinger (1955), who derive it from the pseudo-conservation of the energy- 
momentum tensor T^. 

In this Appendix, for the sake of simplicity, we omit the subscript ^ to 

5 = and 9 = — — — r, (A3) 

Pde a p DE {l + w) 

as well as to the energy density p and to the pressure p of DE. 

The proof is simpler if we start from eq. (A2) and obtain eq. (Al). Let us first 
multiply both sides of eq. (A2) by (1 + w)p and replace 5 by 

^-^ = ^ + 3*(l + «,)^. (A4) 
p p p p a p 

In this way, we obtain 

p(l + w)(e + h = -5p - 3-(l + c 2 s )5p . (A5) 

Owing to the definition of 9 (eq. A3) and thanks to the relation (1 + w)p — 
p+p = (fi 2 /a 2 , the l.h.s. of eq. (A5) becomes (k 2 /a 2 )(j)5(f)+(<f) 2 /a 2 )(h/2), while, 
at the r.h.s., we may replace 5p and 5p with the expressions obtainable from 
the definition of p = 2 /2a 2 + V. Both sides of the resulting equation must 
then be divided by 0/a 2 , so obtaining 



h 2 5(f) + (f)- = -5" ( f)-(i + a 2 ^-2*)5(f)-V"a 2 5(f)-3^(l + c 2 s )(^ + V > 5(l))i . 
2 0a a a 2 (f) 

(A6) 

Here we ought to use the equation of motion of the unperturbed field 0, which 
can be easily reset to yield that 0/0 = —2a/a — a 2 V'/<j). Rearranging the 
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terms, we then obtain: 

8<t> + 2-6'<j> + (k 2 + V"a 2 ) 6</> + £ = 3-[(l - c 2 s )5(f> - (1 + c 2 s )a 2 ^5(f>] . (A7) 
a, 2 a (J) 

Eq. (A7) coincides with the equation of motion (Al), provided that the r.h.s. 
vanishes. This follows the very definition of sound velocity 

2 5p 4>54>/a 2 - V'5(t> ,, 8 , 

from which we easily obtain that (1 + c 2 )/(l — c 2 ) = <f)5<f)/a 2 V'5(f). This last 
relation clearly yields the vanishing of the r.h.s. of eq. (A7). 
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Fig. 1. Optical depth to CMB photons, in different models with fi m = 0.3, h = 0.7 
and n b h 2 = 0.022, as a function of the reionization redshift z r {. The two horizontal 
dotted lines show the reionization redshift consistent with r = 0.08 or 0.16, in the 
case of a single, fast and complete reionization event. 




Fig. 2. Angular spectra of RP models with r = 0.08. Continuous (dashed, dotted) 
lines correspond to A = 10 3 (10 6 , 10 9 )GeV. The absolute values of Q x is plotted; 
changes of sign are indicated by cusps; notice cusps at low /. 
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Fig. 3. As Fig. 2, but for r = 0.16. 
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Fig. 4. As Fig. 2, but for SUGRA potentials. Also for SUGRA potentials, the main 
dependence on the energy scale A concerns the spectrum Cf- at quite low /. This 
dependence is however much milder than for RP potentials. 
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Fig. 5. As Fig. 4, but for r = 0.16. 
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Fig. 6. The likelihood of various energy scales A is determined from artificial data, 
created for a model with given A and for an instrument measuring polarization 
with given a^ ix (indicated in the frames). In these plots A = 10 3 GeV, continuous 
(dashed) lines refer to models with r = 0.16 (0.08). The curves give the likelihood 
obtained averaging among 1000 model realizations. Even with a^ ix ~ 0.1 /uK, there 
is a significant likelihood for A < 1 GeV. See Figs. 11 to relate such likelihood to the 
fraction of realizations for which the model cannot be distinguished from ACDM. 
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Fig. 7. The same as Fig. 6, but for A = 10 6 GeV. Such model is only in marginal 
agreement with data on SNIa. 
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Fig. 8. In this (and the following) plot(s) we describe the effects of cosmic variance, 
showing how 1-a lower limits on A depend on the model realization. The continuous 
(dashed) curve refer to a^ ix = 0.1 (1) /iK. In the frames, the values of the energy 
scale A and of the opacity r are indicated. Notice that, even in the worst case, 
there is a substantial fraction of realizations which allow to set a limit A> 1 GeV, 
so distinguishing a RP model from ACDM. For A = 10 6 GeV and r = 0.16, the 
fraction of realizations indistinguishable from ACDM, with the reasonable noise 
level a^ ix = 1 /xK, is ~ 2 %. 



24 




Fig. 9. As Fig. 8, at the 2-a confidence level. For A = 10 6 GeV and r = 0.16, the 
fraction of realizations indistinguishable from ACDM, with a^ ix = 1 /xK, is ~ 16%. 
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Fig. 10. This figure is analogous to Fig. 6 and concerns the SUGRA model for which 
A is best determined by artificial data. The pixel noise is 10 times smaller than the 
most optimistic expected noise for SPOrt. This plot is obtained by averaging over 
1000 realizations. 
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Fig. 11. 1-cj lower limit distribution for the SUGRA model for which A is best 
determined by artificial data. 



1000 



100 - 



Q 

w 

■A 

I 




Fig. 12. The ratio —h/29DE is plotted with a thick (thin) line for RP models with 
A= 10 3 (10 6 )GeV 
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Fig. 13. Likelihood distribution on £Ide = when RP artificial data are fitted 
with ACDM models. RP model parameters are shown in the frames. In all cases the 
true model has £Ide = 0.7. 




Fig. 14. Ci spectrum of the ACDM model best fitting RP model data with 
A = 10 6 GeV. See text for more details. 
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Fig. 15. Likelihood distribution on Odb, when SUGRA artificial data are fitted with 
ACDM models. In this case, the model parameters, shown in the frames, are well 
approached. 
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Fig. 16. C, spectrum of the ACDM model best fitting SUGRA model data with 
A = 10 3 GeV. See text for more details. In this case, thanks to the lack of negative 
TE correlation, the fitting procedure has no serious misleading effect. 
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